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We study static, spherically symmetric mixed configurations with a nontrivial (wormhole) space- 
time topology provided by the presence of two interacting ghost scalar fields. Wormhole is assumed 
to be filled by a perfect relativistic neutron fluid modeled by a polytropic equation of state. For 
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shown that the maximum of the fluid density is always shifted from the center, and the resulting 
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I. INTRODUCTION 

In recent years interest in obtaining wormhole-type solutions has increased appreciably, primarily because of the 
discovery of the accelerated expansion of the present Universe. At the present time it is widely believed that such an 
acceleration is due to the presence of some special form of matter - dark energy (DE) [l[. One of the key features 
of the latter is its ability to violate various energy conditions. In the most extreme case, DE is modeled by matter 
violating the null energy condition, when the effective pressure p of matter filling the Universe is negative and larger 
(modulus) than its energy density e, i.e., p < —e. Such matter is called exotic. 

In modeling the accelerated expansion of the Universe, it is usually assumed that DE is distributed homogeneously 
and isotropically on the largest scales. This, however, does not exclude a possibility that DE might cluster on relatively 
small scales comparable to sizes of galaxies or even of separate stars. In the latter case, the literature in the field 
offers compact objects composed of DE (the so-called dark energy stars), which are modeled by some form of matter 
possessing the properties of DE 0. Such objects may have both a trivial and a nontrivial topology of spacetime. In 
the latter case, to construct localized solutions, one needs to invoke some form of exotic matter which permits compact 
configurations to have a nontrivial (wormholelike) topology. In the simplest case one can employ the so-called ghost 
(or phantom) scalar fields, which may be massless Q or possess a potential energy Q (other examples of phantom 
field wormholes can be found, e.g., in Refs. Q, and for a general overview on the subject of Lorentzian wormholes, 
see the book Q). 

Another possibility is to consider the case where exotic matter is clustered in the presence of ordinary (baryonic) 
matter. Such mixed compact configurations, with a nontrivial topology provided by using ghost scalar fields, have 
been studied in Refs. [7H10|. In that case, a wormhole is assumed to be filled by ordinary (neutron) matter, and 
the resulting systems possess properties both of wormholes and of ordinary stars: on the one hand, the topology of 
spacetime is nontrivial, and on the other hand, there is a visible surface created by a neutron fluid. Due to the presence 
of the nontrivial topology, the physical properties (masses, sizes, distributions of matter) of such configurations may 
differ considerably from those of systems with a trivial topology (for example, ordinary neutron stars). 

The possibility to construct objects which acquire new properties in the presence of a nontrivial topology motivates 
one to study mixed systems supported by other sources of exotic matter. In the simplest case one can modify the 
systems of Refs. [7 h 10I| “one scalar field plus ordinary matter” by adding to them one more scalar field interacting 
with the first one. Systems with two scalar fields are themselves well known in quantum field theory, where they 
are used to obtain solitary wave solutions [HI- When a gravitational field is present, such systems have also been 
repeatedly considered in the cosmological and astrophysical contexts (see, for example, in Refs, [l^, as well as at the 
beginning of Sec. im. The presence of a second scalar field allows one to get new interesting solutions, not shared 
by one-field systems (for example, the cosmological quintom models containing usual and ghost scalar fields and 
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possessing the properties both of quintessence and of phantom models Q). We will show below that the inclusion in 
our mixed system of the second field will allow us to obtain configurations that possess new physical properties and 
differ substantially from one-field systems. 

Consistent with this, in the present paper we consider mixed astrophysical systems consisting of a wormhole 
supported by two interacting ghost scalar fields and filled by ordinary matter. Our goal is to clarify the question 
of how the presence of such a wormhole changes the distribution of the ordinary matter and influences the physical 
characteristics of the resulting objects. 

The paper is organized as follows. In Sec.|n]we present the statement of the problem and derive the corresponding 
general-relativistic equations for the mixed systems under consideration. In Sec. IIIII we solve these equations nu¬ 
merically for different amount of ordinary matter and compare the systems under consideration with configurations 
studied earlier. Finally, in Sec. 113 we summarize the results obtained. 


II. STATEMENT OF THE PROBLEM 


As pointed out in the Introduction, here we study a mixed gravitating system consisting of two interacting ghost 
scalar fields (f, x a perfect fluid. The Lagrangian for this system can be presented in the form 


L = -- 


IGttC? 


zR- 


+ ^d^xd^X - V(¥’, x) 


+ Rri 


( 1 ) 


with the curvature scalar R and Newton’s constant G. Here is the Lagrangian of the perfect isotropic fluid (where 
isotropic means that the radial and the tangential pressure of the fluid agree), which has the form Lm = p. 

We choose the scalar field potential energy V {p, x) as 

yi.^,x) = -mlf + ^{x^ -mlf + ip^x^ -Vo- ( 2 ) 

Here mi and m 2 are the masses of the scalar fields, Ai,A 2 are the coupling constants, and Vq is some constant 
whose value is chosen from the statement of the problem. (Note that one of these free parameters can always be 
eliminated by the corresponding rescaling.) Using this potential, we obtained earlier a number of solutions which can 
be employed both in describing astrophysical objects and when considering cosmological problems. Namely, we have 
shown that: (a) for the four-dimensional case, there exist regular spherically and cylindrically symmetric solutions 
[mil, and also cosmological solutions both for usual fields (i.e., for the fields having the usual sign in front of the 
kinetic energy term) and for ghost scalar fields [H : (b) for the higher-dimensional cases, there exist thick brane 
solutions supported by usual and ghost scalar fields [l7| . 

An important feature of this potential is the presence of two local minima at y = 0, (^ = ± mi to which the solutions 
tend asymptotically at spatial infinity. These local minima correspond to two vacua. For wormhole-type systems, the 
solution starts in one of these vacua as r —>■ — 00 and returns to it back as r —>■ -l-oo [l^ (r is the spatial coordinate). 
Such solutions are called nontopological m , in contrast to solitonlike solutions with one scalar field when the solution 
starts in one of the vacua and goes to the other. The latter solutions occurring only in the presence of two or more 
vacua are called topological. In the present paper we will discuss only the nontopological solutions. 


A. Field equations 

To describe the mixed equilibrium system under consideration, let us employ the polar Gaussian coordinates, in 
which the metric has the form 


ds^=e‘'{dx^) —dr^ — R^dXl^, 


(3) 


where v and R are functions of r only, and = ct is the time coordinate. 

The total energy-momentum tensor can be obtained from the Lagrangian 0 in the following form: 

Tt = ie + p)u^u’^ - S’^p - d.ipd^ip - dtxd^x - - 7 ;d^,xd^X - V(v>, x) 


(4) 
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where e and p are the energy density and the pressure of the fluid, and u* is the four-velocity. By using ([4]), the (q), 
(}), and (2) components of the Einstein equations with metric ([3]) take the form 


[2^" + 

fR'y 


KRJ _ 


i?2 


SttG. 


.4 ^0 


rpO _ 

-in — 


SttG 


e - 2 + -yi^,x) 


~~R 


1 SttG 1 SttG 


-P+ 2 + X'^) -V(‘P,X) 


R" IR', 1 ,, 


= - 


SttG , SttG 
—= —- 


P+2 + X'^) +1^(+>X) 


(5) 

( 6 ) 
(7) 


where the prime denotes differentiation with respect to r. 

The equations for the scalar fields p, x result from the Lagrangian (ED in the general form as 


_^ 

^/^ dx>^ 


9 


flV 


d{p,x) 

dx^ 


In the metric (O, these equations yield 


dV 

xY 


( 8 ) 




X 


1 / 

-u' -f 2— 
2 R 

1 / 

2*^ 


if = 


x' = 


dV 
dp ’ 

dx ' 


(9) 

( 10 ) 


Finally, the hydrostatic equation for the fluid can be obtained from the law of the conservation of energy and 
momentum, = 0. The i = 1 component of this equation gives 


dTl 1/1 n\ / dd! 

1F+2 + -T»).' + 2- 


r;-l(r| + T|) 


= 0 . 


( 11 ) 


Taking the components of the energy-momentum tensor from Eqs. dSD-ED, and also taking into account that Tg = T2 , 
from EID we have 


dp 

dr 


1 , ^dv 

-2'" + '’+' 


( 12 ) 


The mixed system under investigation can be regarded as follows. Let initially we have a wormhole geometry 
provided by two interacting ghost scalar fields. In the absence of ordinary matter (fluid), such systems have been 
considered in M- By adding a fluid to the system, here we consider the influence that such a wormhole has on the 
spatial distribution of the fluid and on the physical characteristics (masses, sizes) of the resulting mixed systems. 

As the fluid one can choose any type of matter used in modeling astrophysical objects (for instance, stars). One of 
the simplest ways to describe star’s matter is to employ a polytropic fluid. The latter is applied both in describing 
nonrelativistic objects (Newtonian stars, see, e.g., [T^) and when considering relativistic systems, including those 
with a relativistic neutron fluid in a strong gravitational field. In doing so, a polytropic equation of state (EOS) can 
adequately represent more realistic EOSs used in modeling neutron stars (see, e.g., |19|b 
For our purpose, we choose the following relativistic EOS: 


p = Kp 


1+1/n 
b > 


e = pbC^ + np, 


(13) 


with the constant K = kc^rribY f^^e polytropic index n = 1/(7 — 1), and pb = Ubnib denotes the rest-mass 

density of the neutron fluid. Here rib is the baryon number density, is a characteristic value of rif,, mb is the 

baryon mass, and k and 7 are parameters whose values depend on the properties of the neutron matter. In particular, 
here we choose mb = 1.66 x 10“^^ g, = 0.1 fm“^, k = 0.1, and 7 = 2 [2^. 

With the EOS in the form of (IT3l) . one can integrate Eq. EID- +0 ‘io f^is, it is convenient to introduce the new 
variable 9 (l^ . 


Pb = PbcO^ , (14) 

where pbc is some characteristic density of the neutron fluid. For the mixed configurations with an isotropic fluid 
considered by us earlier in Refs. [3, 0 Pbc corresponds to the central (maximum) density at the wormhole throat. 
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However, for systems with an anisotropic fluid a situation is possible where the maximum density of the neutron 
fluid is shifted away from the center |10l |. Then no longer plays the role of maximum density, but it is just some 
characteristic value for the configuration under consideration. As we will see below, here we deal with such a situation 
even in the case of an isotropic fluid. 

Making use of expression o, we have from Eq. (USD 

2a{n +1)^ =-[1 + a{n + 1)9]^, (15) 

dr dr 

where a = jc^ is the relativistic parameter [^ . Integrating this equation, one can find 

^ 

where e‘'‘= is the value of e'^ at the center. The integration constant Vc is fixed by requiring that the spacetime is asymp¬ 
totically flat, i.e., = 1 at infinity. In turn, the arbitrary constant 0c, contained in the boundary conditions (I24L 

corresponds to the central value of the function 9. 

For the numerical calculations, it is convenient to rewrite the obtained equations in terms of dimensionless variables. 
Since below we will seek solutions with zero central values of the scalar-field derivatives, <p^(0), = 0, ^ nonzero 

central value of the potential energy, l/((p(0), x(0)) ^ 0, we normalize the characteristic size of the system L with 
respect to this value. Namely, let us introduce dimensionless variables: 

r R 

^=7> where L = r—^ =■ (17) 

L L c2 v'8^G|H(^(0),x(0))| 

In such variables Eqs. ©-dzi), dH), (uni) take the form 


/E' 

2 - 1 - — 

E V S 

E' /E' 


+ ^ = H(1 + an9)9- - i + x'^) - V, 


EVE 


+ ^' ) + ^ = -Bct 9-^^ + \ {(p'^ + X'") - F, 


¥ + +r'" = ^ 

(p" -I- ‘p’ = ^ [2x^ + Ai {(p^ — , 

r + + 2 ^) x' = X [2^^ + A 2 (x^ - pI)] , 


where the dimensionless potential energy is 


plf + ^(x^ - plf + p^x^ - %. 


(18) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 

(23) 


Here B = (pf,cC^)/|17((p(0), x(0))| is the dimensionless ratio of the characteristic fluid energy density to that of 
the scalar field at the center; the dimensionless coupling constants Ai _2 = ^c'^/87rG-\/|H((p(0),x(0))|^ Ai^ 2 , and 

l^o = l^)/in¥^(0),x(0))|. 

Note that, using L from 113, the expression for B can be recast in the following form: B = ^ttG pbc{L/cY. This 
expression will be used below in performing numerical calculations. In the case of H = 0 we deal with a system 
consisting of a pure scalar field configuration with no ordinary matter [T3 | . 


B. Boundary conditions 

Here we consider neutron-star-plus-wormhole configurations that are asymptotically flat and symmetric under 
^ —>■ —The metric function E(^) may be considered as a circumferential radial coordinate. Asymptotic flatness 
requires that E(^) —>• |^| for large |^|. Because of the assumed symmetry, the center of the configurations at ^ = 0 
should correspond to an extremum of E(^), i.e., E'(0) = 0. If E(^) has a minimum at ^ = 0, then ^ = 0 corresponds 
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to the throat of the wormhole. If, on the other hand, has a local maximum at ^ = 0, then ^ = 0 corresponds to 
an equator. In that case, the wormhole will have a double throat surrounding a belly (for examples of double-throat 
systems, see, e.g., Refs. [fli. 1^1. 

Consistent with this, Eqs. (ra-®, together with the expression (fTOl) . will be solved for given parameters of the 
fluid cr, n, and B, subject to the boundary conditions at the center of the configuration ^ = 0, 

0(0) = E(0) = Ee, S'(0)=0, z^(0) = ^c, i^'(0)=0, m = <p'(0) = 0, x(0) = Xc x'(0) = 0, 

(24) 

where 0ciSc,(^ci and Xc are some constants. 

To find the location of the throat, expand the metric function E in the neighborhood of the center as 

E«E, + l/2E2e". 

Then, using Eqs. (fTSl) and (IT^ . we hnd the relations 

Ec = — , E2 = —[1 -|- (n -|- l)a0c] ■ (25) 

V-Sa0c”+'-E(^c,Xc) 

It is seen from these expressions that (i) to obtain physically reasonable solutions, one needs to choose the free 
parameters of the system in such a way that the radicand in the first formula be positive, and (ii) at 0^ > 0 (as it 
should be for ordinary matter) E2 is negative, corresponding to a solution with an equator surrounded by a double 
throat. 


III. NUMERICAL RESULTS 

In Ref. [13 we studied wormhole systems supported by ghost scalar fields with the potential ([2]). Here we modify that 
system by adding to it ordinary matter in the form of a polytropic fluid. Our aim is to study such mixed configurations, 
in particular, with respect to the effect that the presence of the nontrivial topology has on the distribution of the 
neutron matter along the radius of the system. 

We will solve the system of equations dUD-® together with (ITSl) numerically using the boundary conditions (IMl) 
and (1251) . Since a characteristic property of a gravitating polytropic fluid is the presence of an edge where its pressure 
(density) goes to zero, the configurations under consideration can be subdivided into two regions: (i) the internal one, 
where both the scalar fields and the fluid are present; (ii) the external one, where only the scalar fields are present. 
Correspondingly, the solutions in the external region are obtained by using Eqs. (HU-dlii), in which 0 is set to zero. 

The internal solutions must be matched with the external ones at the boundary of the fluid, C = Cfej by equating 
the corresponding values of the functions <f,x^ ^ their derivatives. The boundary of the fluid is defined 

by p{^b) = 0. Knowledge of the asymptotic solutions in turn allows one to determine the value of the integration 
constant Vc a-t the center, proceeding from the requirement of asymptotic flatness of the external solutions. 

As in the case without ordinary matter studied in Ref. , the system of equations (rT^ - (l22ll has regular solutions 
(satisfying the necessary boundary conditions) only for certain values of the coupling constants Ai,A2 and of the 
masses of the scalar fields pi, fi2- As a result, the problem reduces to a search for eigenvalues of the parameters pi, p2 
and for the corresponding eigenfunctions i/, E,:^, and x of the nonlinear system of differential equations (11811 - (I22[) . 
We will seek the specified eigenvalues by u sing the shooting method. A step-by-step description of the procedure for 
finding solutions can be found in Refs. |14l. 11^. 

Proceeding in this way, we have obtained the results shown in Figs. [TJ131 The parameters of the system have been 
chosen so that the masses of the mixed configurations under consideration be of the order of 1-2 solar masses, and 
their sizes ^ 20 km. These characteristics are typical for neutron stars . 

The results were obtained by choosing the following fixed values of the parameters: the characteristic fluid density 
Pbc = 10^'^g cm“^ and its size L = 10 km, the coupling constants Ai = 0.1, A2 = 1, the central values of the scalar 
fields (pc = 5,Xc = l-l- This allowed to get configurations with the required characteristics, as shown in Table HI 

As seen in Fig. ^ (p ^ and x 0 at large distances, i.e., the solutions approach asymptotically the local 
minimum of the potential (E^ . The arbitrary constant 1^ is chosen so that the value of the potential (1251) be equal 
to zero at the local minimum, i.e., we take Vq = A2/r|/4. Such a choice ensures a zero value for the scalar field energy 
density as ^ ±00, as one can see from Fig. [TJ 

The asymptotic behavior of the scalar fields can be found in the following form: 


ip = p.i-6>p, X = Sx, 


( 26 ) 
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FIG. 1: The fluid energy density B{1 + Gn9)9^ (left panel) and the total energy density Tq (right panel) from the right-hand 
side of Eq. mu (in units of |E((p(0),x(0))|) are shown as functions of the relative radius ^/^b- Since the solutions with the 
wormhole are symmetric with respect to ^ = 0, the graphs are shown only for ^ > 0. In the left panel, the parameter 9c runs 
the values 0, 0.1, 0.4, 0.7,1.0,1.5, from bottom to top. In the right panel, all the remaining curves with 0 < < 1.5 lie between 

the presented curves. The thin vertical line corresponds to the boundary of the fluid. 


TABLE I: Characteristics of the mixed configurations at fixed pbc, L, Ar, A 2 , ipc, Xc (see in the text). Here the central pc = Pbc9c 
and the maximum pmax densities of the neutron fluid (both in units of 10^‘^g cm“®), the total mass Mtot, the proper mass 
Mprop of the neutron matter (both in solar mass units Mq), the radius of the neutron fluid R (in kilometers), and the masses 
of the scalar fields pi,p 2 are shown. For the neutron star with the same values of the polytropic parameters and pbc, we have 
Mtot = 1.876M0 and R = 32.812 km. 


Pc 

Pmax 

Mtot 

-^^prop 

R 

pi 

P2 

0.0 

0.838 

1.087 

0.526 

18.555 

5.97803 

7.67266 

0.1 

0.932 

1.173 

0.623 

19.201 

5.97824 

7.67252 

0.4 

1.213 

1.419 

0.902 

20.437 

5.97884 

7.67213 

0.7 

1.489 

1.631 

1.147 

20.980 

5.97940 

7.67175 

1.0 

1.761 

1.805 

1.351 

21.132 

5.97994 

7.67140 

1.5 

2.206 

2.022 

1.613 

20.930 

5.98078 

7.67084 


where 6(p, ^ 1 behave as 


exp(-^/2A^n exp - A 2 Mi) 

^^^^( 27 ) 

where C^p, are integration constants. 

The total masses of the configurations shown in the Table were calculated as follows. Consider a volume enclosed 
by a sphere with circumferential radius Rc, corresponding to the center of the configuration, and another sphere with 
circumferential radius R> Rc- Then, using the spherically symmetric metric ([3]), the mass m(r) associated with this 
volume can be defined as follows: 


= —Rc 


Att 


T^R^dR. 


( 28 ) 
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FIG. 2: The metric functions gu = e.’' (left panel) and E (right panel) are shown as functions of the relative radius C/Cfc- 
the left panel, the thin vertical line corresponds to the boundary of the fluid; the dashed line refers to the neutron star; for the 
mixed configurations, the parameter 6c runs the values 0,0.4, 0.7,1.0,1.5, from top to bottom. Asymptotically, as ^ ±oo, 
the spacetime is fiat with E ^ |^| and ^ 1 from below. 



FIG. 3: The typical behavior of the scalar fields. Asymptotically, as f —^ ±oo, the field (3 tends to ui and y goes to 0 (see 
TableUand Eqs. ([MI-llZl]. 


For single-throat configurations, the circumferential radius Rc corresponds to the radius of the wormhole throat 
dehned by Rc = min{i?(r)}. However, for the systems considered here, a single-throat configuration will exist only 
when Oc = 0. When 6c ^ 0, as mentioned above, Rc corresponds to an equator, while the two throats are located 
symmetrically away from the center (see in the right panel of Fig. El). 









In the dimensionless variables of Eqs. (II3 and (fl3) the expression (|28l) takes the form 

>dl 


m(0 = M* + /^ \b {1 + ane) 9 - - i + x'^) - v] 


(29) 


where the coefficient M* in front of the curly brackets has the dimension of mass 

B 

SirG^Pbc' 

Note that formally the total mass Mtot is then obtained by taking the upper limit of the integral to infinity, since 
the energy density of the scalar fields becomes equal to zero only asymptotically, as E —>■ oo. However, since the scalar 
fields decay exponentially, see Eq. (1271), the integral converges rapidly to its asymptotic value even at small ^ ^ 0(1)- 
Note also that in evaluating the above integral it is necessary to perform the calculations in the internal and external 
regions separately. 

One more interesting characteristic of the configurations under consideration is the proper mass Mp^^p of the neutron 
matter: 

fi"b rG 

Mprop = rribN = An j pbR^dr = AnpbcL^ / d"E^d^. (30) 

Jo Jo 

Mprop is equal to the mass which the baryons of the system would possess altogether, if they were dispersed throughout 
a volume so large that all types of interactions between them could be neglected. Evaluating the expression (1501) for 
the proper mass, we can find the number of neutrons N in the system. 

From the results obtained, it is possible to delineate the following features of the systems under consideration. 

(1) For the values of the parameters used here, it is seen from the Table that while the central value of the neutron 
fluid density increases, the masses of the scalar fields remain practically unchanged. This corresponds to the fact 
that the spatial distribution of the scalar fields changes also very slightly for different amount of the neutron 
matter. In this sense, the scalar fields can be regarded as the background ones. But the total mass of the 
configuration varies substantially that is obviously caused by a change in the number of neutrons N in the 
system (see the Table). 

(2) The presence of the ordinary matter results in the following changes in the geometry of the system: (a) as the 
amount of the neutron matter increases, the redshift function gu = increasingly differs from 1 (see Fig. Ill), 
i.e., the system becomes more and more relativistic; (b) at 0c ^ 0 two weakly marked throats are necessarily 
present in the system (see in the right panel of Fig. [2|), in contrast to systems supported by scalar fields only 
where just one throat is present. 

(3) One more interesting feature associated with the presence of the two scalar fields is a unique shape of the 

neutron matter distribution along the radius. The literature in the field offers different systems with an isotropic 
fluid having a maximum density at the center, like that illustrated by the curve labeled by “neutron star” in 
Fig. □ This applies both to systems with a trivial topology of smcetime (ordinary stars) and to the mixed 
systems of the “wormhole plus ordinary matter” type of Refs. When a system contains an anisotropic 

fluid, its maximum density may already be located somewhere between the center and the edge of the system 
(see Refs. [2^ concerning ordinary stars and Ref. [13 where a mixed system is discussed). For the mixed 
configurations considered here such a shift of the maximum fluid density takes place even for isotropic matter, 
and the central density of the neutron matter can vary from 0 to the values of the order of nuclear density (see 
the Table and Fig. [T]). 

Of course, the obtained quantitative characteristics (masses, sizes, shapes of the neutron matter distribution) 
depend essentially on the specific values of scalar field parameters (coupling constants, boundary conditions) and 
of ordinary matter characteristic densities and sizes. But one might expect that the qualitative behavior of the 
solutions will remain the same even with a different, physically reasonable choice of the parameters of the system 
under consideration. 



IV. CONCLUSION 

We have studied equilibrium mixed configurations consisting of a wormhole supported by two interacting ghost 
scalar fields and threaded by a relativistic polytropic neutron fluid. The scalar fields play the role of exotic matter 
providing a nontrivial topology of spacetime in the system. 
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Our goal was to investigate the influence that such a nontrivial topology has on the distribution of the neutron 
matter and on physical characteristics of resulting objects. In doing so, we compared the obtained mixed configurations 
with a neutron star modeled by matter with the same EOS. Having fixed the polytropic parameters and coupling 
constants of the scalar fields, we kept track of changes in the distribution of the neutron matter depending on the 
central value of the neutron matter rest-mass density. 

The main results can be summarized as follows: 

(i) There exist static regular asymptotically flat solutions describing mixed neutron-star-plus-wormhole systems 
in which the neutron matter is concentrated in a finite-size region. Such configurations may be regarded as 
consisting of a neutron star having either a single throat at the center (at zero central density of the neutron 
matter) or an equator surrounded by a double throat (see Fig. [2|). 

(ii) Due to the presence of two scalar fields, the neutron fluid always has a maximum density somewhere between 
the center of the configuration and the edge of the fluid. In this case the central density of the neutron matter 
can vary from 0 to the values of the order of nuclear density (see Table Hand Fig. [T|). Such a behavior, which 
is usually characteristic of anisotropic fluids, is attained here even in the case of isotropic fluid. This is because 
of the characteristic behavior of the redshift function gtt = which, in systems with two scalar fields, has a 
minimum shifted away from the center (see Fig. [2 and also the pure scalar field systems of Ref. H). 

Altogether, as regards the form of the neutron matter distribution, the mixed systems “two scalar fields plus isotropic 
matter” discussed here are similar to the mixed systems “one scalar field plus anisotropic matter” of Ref. m , with 
the main difference related to the behavior of the redshift function. On the other hand, the mixed two-field systems 
considered here differ in principle from the mixed one-field systems of Refs. [3, Q, for which the maximum of the 
neutron fluid density is always located at the center of a configuration. 

In conclusion, we may say a little regarding the question of stability of the considered mixed systems. The results 
available in the literature indicate that for pure field configurations (with no ordi nary matter) wormhole solutions 
are linearly [25l - [^ and nonlinearly unstable (see, however, the recent work (30| where a method of obtaining 

general-relativistic stable wormhole solutions is presented). As shown in our previous studies of mixed “wormhole 
plus star” systems with one scalar field [I®, the presence of ordinary matter did not result in stabilization of the 
configuration. One might expect that the inclusion in the system of one more ghost scalar field will not change the 
stability situation. Nevertheless, this issue requires special study, and we plan to perform such a stability analysis 
later on. 
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